Temperature and doping dependence of high-energy kink in cuprates 
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It is shown that spectral functions within the extended t-J model, evaluated using the finite-temperature 
diagonalization of small clusters, exhibit the high-energy kink in single-particle dispersion consistent with recent 
angle-resolved photoemission results on hole-doped cuprates. The kink and waterfall-like features persist up to 
large doping and to temperatures beyond J hence the origin can be generally attributed to strong correlations 
and incoherent hole propagation at large binding energies. In contrast, our analysis predicts that electron-doped 
cuprates do not exhibit these phenomena in photoemission. 

PACS numbers: 71.27,+a, 75.20.-g, 74.72.-h 



The anomalous properties of quasiparticles (QP) in 
cuprates are most directly probed by the angle-resolved pho- 
toemission spectroscopy (ARPES) [1]. Besides the Fermi sur- 
face development, pseudogap features and low-energy kink 
|0]] considerable attention has been recently devoted to the 
high-energy kink (HEK) observed by ARPES quite univer- 
sally in hole-doped cuprates lHUHBSSSSEl. The 
anomaly appears in the QP dispersion along the zone diag- 
onal (0,0) - (7r/2,7r/2) as a kink at binding energies typi- 
cally Ei ~ 0.4 eV followed by a fast drop - 'waterfall' to 
the high-energy scale E2 ~ 1 eV being well pronounced 
around the zone center k ~ (0,0). The HEK seems to ex- 
ist in a broad range of hole-doped cuprates: from undoped 
Jl H, underdoped SHU, optimally doped fl 11 SB, 
overdoped JH H 0, HI] to highly overdoped regime yl 9]. 
Some ARPES spectra H 01 give indications for the coexis- 
tence of two branches in the part of the Brillouin zone: one 
representing the renormalized QP band reaching E ~ E\ and 
the second being the remnant of the high-energy unrenormal- 
ized band at E ~ E2. Analogous evidence for the HEK in 
electron-doped cuprates [7j,ll0|] is less conclusive. 

Theoretically, the origin of the HEK is presently lively de- 
bated. Since in contrast to low-energy kink (H the energy 
Ei is too high to be attributed to phonons, several aspects of 
strong correlations are given as a possible explanation. The 
similarity to spectral functions of one-dimensional (ID) chain 
cuprate SrCu02 ifTlll with pronounced two component spec- 
tra, i.e., spinon and holon branches, seems to support the long- 
sought spinon-holon scenario also for two-dimensional (2D) 
cuprates [4] . On the other hand, alternative explanations with 
string excitations of a QP in an antiferromagnet (AFM) 111 211 . 
split QP band within the slave-boson theory 11311 . and the 
vicinity to a quantum critical point H 1 411 are not restricted to 
ID. Recent numerical calculations within the Hubbard model 
support the existence of the HEK in prototype models of cor- 
related electrons, both for the undoped system lfl5ll as well as 
in the large-doping regime 01611 . where the origin of HEK is 
attributed to high-energy spin correlations 1 lrjl 171. 

In the following we present finite-temperature numerical re- 
sults within the prototype t-J model of strongly correlated 



electrons in cuprates. They reveal the existence of the HEK 
in a broad range of hole concentration Ch and temperature T 
in the uj < sector of spectral functions A(k, oj), correspond- 
ing to ARPES in hole-doped cuprates. Well pronounced at 
intermediate and large doping as the waterfall-like dispersion, 
the HEK develops at lower doping and T < J into two partly 
coexisting branches, the renormalized QP band and a broad 
bottom band. An important fact for the interpretation is the 
observed persistence of the HEK up to high T ~ t > J which 
gives strong support to the scenario that the HEK and waterfall 
are quite universal signatures of strong correlations and only 
indirectly connected to low-T phenomena as the longer-range 
AFM and superconductivity in these materials. Also, strong 
asymmetry in li leads to the conclusion that analogous phe- 
nomena in electron-doped cuprates should be absent within 
ARPES spectra. 

We study the single-particle excitations within the extended 
t-J model 
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where c| are projected fermionic operators not allowing for 
the double occupancy of sites. As relevant for cuprates we 
consider the model on a square lattice and include besides the 
nearest-neighbor = t also the second-neighbor tij = t' 
and the third-neighbor hopping = t". We present in the 
following results for t' = -0.25 t, t" = 0.12 t, J = 0.4 t 
iflii Il9tl chosen to reproduce well properties of hole-doped 
cuprates, e.g., the measured Fermi surface. 

We calculate the spectral function A(k, u>) using the 
usual T — exact diagonalization method and the finite- 
temperature Lanczos method (FTLM) for T > [20]. Sys- 
tems considered here are tilted square lattices of N = 18, 20 
sites with finite concentration of holes ct = Nh/N doped 
into the reference undoped AFM insulator. Since fixed bound- 
ary conditions on small systems allow only a discrete set of 
wavevectors ki; , / = 1, N we employ twisted boundary con- 
ditions to scan the whole Brillouin zone 1 19], k = k; + 9 by 



introducing hopping elements £j 



tij — tij exp(i6 



Eq. fl}. For details of the application of the FTLM to spectral 
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functions we refer to Ref . J2lll . Besides the evident possibil- 
ity of obtaining T > results the FTLM allows for a reliable 
evaluation of the self energy S(k, uj) which is essential for the 
interpretation of observed phenomena. 

First we present results for A(k, uj) at hole-doping ct = 0.1 
corresponding to underdoped regime calculated on a system 
of N = 20 sites. In Fig. 1 we present the weight map of 
A(k, uj) along the diagonal and the edge directions within the 
first Brillouin zone and its evolution with increasing T. The 
T = result in Fig. 1 a is obtained by the ground state Lanczos 
procedure 01911 . while in Fig. lb,c,d FTLM results are shown 
for increasing T/t = 0.2, 0.4, 0.75. 
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Figure 1: (Color online) Weight map of A(k, ui) vs. k along sym- 
metry lines in the Brillouin zone for cu = 0.1 and different T/t. 

If one concentrates on the dispersion along the zone di- 
agonal (0, 0) — (71", tt) it is easy to recognize the HEK fea- 
ture at uj ~ — t for all presented T. In fact, a pronounced 
waterfall-like single band dispersion is evident even at very 
high T ~ t > J, where the steep drop appears close to 
k ~ x(ir, tt) with x ~ 0.3. It should be, however, noted that 
T ~ t represents already very high T in this doping regime 
which leads to a substantial shift of the chemical potential so 
that the Fermi surface is tending towards (tt, tt) as evident in 
Fig. Id. 

With lowering T < J the uj = (Fermi surface) cross- 
ing of the dispersion along the zone diagonal approaches 
k ~ (tt/2, 7t/2) as expected for low doping. More relevant 
here, the single dispersion curve evolves into a more complex 
structure: a) the uj > part not accessible by ARPES re- 
veals a well-defined dispersion of weakly damped QP, b) the 
renormalized band with small QP velocity remains well de- 
fined close to the Fermi surface, i.e. at k ~ (tt/2,tt/2), or 
even extends nearly to (0, 0) at low T — > 0, c) less coher- 



ent band-like feature corresponding roughly to the bottom of 
the unrenormalized band is well developed close to the zone 
center, k ~ (0, 0). For T < J both bands coexist at least at 
k - (7r/4,7r/4). 

All observed features are present also in (tt, 0) — (0, 0) di- 
rection as clearly seen in Fig. 1. However, the difference ap- 
pears at higher T where the entire band in this direction is 
positioned well below the chemical potential and does not ex- 
perience the waterfall effect anymore. This is in agreement 
with the evolution of the band towards the usual although still 
renormalized tight-binding dispersion. 
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Figure 2: (Color online) Weight map of A(k, uj) vs. k along symme- 
try lines for fixed T/t = 0.2 and different hole dopings cu- 

Quite analogous behavior can be followed at fixed low T as 
a function of doping. We show results obtained using FTLM 
for systems with N = 18, 20 within a broad doping range 
c h = 0/20 - 4/18. Note that for the undoped system, Ch = 0, 
the absolute position of the chemical potential is not well de- 
fined within the t-J model. In the latter case in Fig. 2a, one 
can again recognize very well pronounced renormalized QP 
band reaching the zone center, while the bottom band is very 
incoherent. Both bands coexist in (ir/4, 7r/4) — (0,0) and 
(7r/4,0) — (0,0) regions. With increasing Ch the renormal- 
ized QP band remains well defined near the Fermi surface at 
k ~ (tt/2, 7t/2), but dissolves towards the zone center. At 
the same time the bottom band starts to move away from the 
renormalized band and stays well pronounced near k ~ (0,0). 
Both bands are connected with the waterfall-like drop of low 
intensity. While the whole effective bandwidth is weakly re- 
duced Aw ~ 6 t < 8 t relative to a tight-binding band, the 
bottom of the band at k = (0, 0) is deeper than expected from 
the tight-binding dispersion in the case of = 3/20,4/18. 
This is consistent with experimental observations yl. 
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The origin of the HEK can be best analyzed and understood 
by expressing the single-particle Green's function correspond- 
ing to A(k, uj) = — ImG(k, oj) in terms of the self energy 

S(k,w), 



G(k,w) 



Ck-S(k,w) 



(2) 



The model, Eq.([T]i, defined with projected fermionic operators 
requires a nonstandard normalization a as well as a nontrivial 
'free' term Ck representing the first frequency moment of the 
A(k,uj) 12211 . Within the paramagnetic metal with (Si) = 
one can express explicitly a = (1 + Cf l )/2 and 



Ck = C-4^^*i7i(k), 



1 



{So -Si), (3) 
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where tj,j = 1,3 represent hopping parameters t, t', t", re- 
spectively, which are renormalized with rj that involve local 
spin correlations (Sq-Sj). The tight-binding band dispersions 
corresponding to ti are then 71 (k) = (cosfc x + cosfc y )/2, 
72 (k) = cosk x cosky and 73 (k) = (cos 2^ + cos2/cj,)/2. 
The above expression, Eq.(f2|, in terms of a and (k leads to 
properly analytically behaved £(k, uj — > ±00) oc 1/oj. 

Following Eq.© we extract £"(k, uj) provided that 
A(k, oj) are smooth enough which is for available systems 
typically the case for T/t > 0.1. In Fig. 3 we present results 
for E"(k, uj) corresponding to spectra in Figs. 1 at Ch = 2/20 
and various T/t, but fixed k = (tt/A, 7t/4) chosen to repre- 
sent the location of the HEK. It should be noted that £" (k, oj) 
is not crucially dependent on k (ignoring here more delicate 
phenomena as the pseudogap [21] ), at least not inside the 
Fermi volume so results in Fig. 3 are representative for all k 
relevant for effective ARPES bands. 

Several characteristic properties of the QP damping recog- 
nized already in previous studies 11201 fell 12211 can be deduced 
from Fig. 3. a) The damping function E"(k, oj) is very asym- 
metric with respect to the Fermi energy oj = 0. For the hole- 
doped case discussed here the damping is large only for oj < 
corresponding to ARPES. b) As one expects in a metal we find 
E"(k, oj = 0) — > (or at least decreasing) at low T — > 0, a 
prerequisite for a well defined Fermi surface, c) Within quite 
a large regime — 2t < oj < we recover at low T well known 
marginal variation — S"(k, ui) oc \oj\ |20[ 23], while only at 
large oj < — 3 t the damping decreases and loses intensity, e) 
Increasing T mainly influences the behavior close to oj ~ 
filling the dip and increasing |S"(k, w ~ 0) |, at the same time 
making S"(k, oj) more featureless. 

Clearly, the strength and the form of E"(k, oj) determines 
the anomalous dispersion oj-^ via the pole location 



1 f ,£"(k,a/) 

- Ck + - / doj = 0. 

n J uj k -oj' 
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The relevant quantity to estimate the influence of E" (k, uj) 
on the dispersion oj^ is the intensity rfc = — j E"(k, oj)duj/ir. 
We notice that at low doping rfc is not strongly dependent ei- 
ther on k, Ch nor on T. In the range of interest corresponding 



Figure 3: (Color online) Damping function — E"(k, ui) correspond- 
ing to Fig. 1 for Ch = 0.1, k = (7r/4, 7r/4) and various T/t = 
0.1-1. 



to Figs. 1-3 we find rfc ~ 3 — 4 t 2 . It should be noted that the 
origin of large 77k is here entirely in strong correlations, i.e., in 
the incoherent motion of a particle (hole) in a spin background 
with singly occupied sites. Such physics can be well captured 
by, e.g., a retraceable path approximation B24I1 where one gets 
rf = At 2 , very close to our numerical results. 

Since £k in Eq.(0) produces only a regular although renor- 
malized tight-binding dispersion, the anomalous effective dis- 
persion emerges from E'(k, oj). Due to large r/^ and a re- 
stricted range —6 t < oj < of appreciable E"(k, ui)\, 
E'(k, uj) leads to a substantial change of the dispersion in this 
oj regime. At low doping and T < J it induces in combination 
with a narrow (k a coexistence of renormalized QP band and 
the bottom band at k < (7r/4, 7t/4). The latter one is quite 
incoherent due to large |E"(k, oj)\ in oj < region. On the 
other hand at uj ~ one has E"(k, uj) — > which allows for 
a well defined renormalized QP band near the Fermi surface. 

The effect of T > is to broaden E"(k, uj) and to in- 
crease QP damping at oj — > 0. Then, E'(k, uj) shows less 
structure and the renormalized QP peak at low oj transforms 
with increasing T into a single effective band. However, due 
to T-independent r/^ the structure of E'(k, uj) remains strong 
enough to keep the waterfall drop up to very high T. Note 
that even for larger T as in Figs. lc,d an effective dispersion 
following Eq.© remains renormalized by rj ~ a although 
the band becomes wider as spin correlations loose intensity 
for T > J. Analogous are phenomena at larger doping except 
that E"(k, uj) generally decreases with c^. 

To illustrate that above features are essential and sufficient 
to reproduce the HEK and the waterfall we compare numeri- 
cal results in Fig. 1 with a simplified model of |E"(k, oj)\ as- 
suming: a) E"(o;) is local, i.e., k independent, b) at T = it 
follows marginal behavior (linear in oj) for — e a < uj < B23I1 . 
c) for larger binding energies —65 < uj < —e a it decreases lin- 
early to zero, d) the effect of T > is to convolute E"(u;, T — 
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0) with usual thermodynamic factor f(u>)[l — /(<*>)] where 
f(u>) is the Fermi-Dirac distribution. 

For results presented in Fig. 4 we fix e a — t, = 6 t, rj — 
2t and vary through r\ = 0.35,0.5 for T = 0, 0.75, respec- 
tively, while T2.3 = a. We present in Fig. 4 the T-dependence 
of A(k, uj < 0) with k along the zone diagonal. It is well 
visible how the two-band structure at low T transforms into a 
rather regular but broad single band with persistent waterfall 
even at very high T ~ t. 



T/t=0.75 




Figure 4: A(k, u) along the zone diagonal calculated from a simpli- 
fied model and various T/t = 0, 0.4 and 0.75. 

In conclusion, we have shown that the prototype model 
as the extended t-J model incorporates the physics of the 
HEK as well as the waterfall as observed in numerous recent 
ARPES studies of hole-doped cuprates. While at low Ch and 
low T the spectra typically reveal a coexistence of a narrower 
renormalized QP band and an incoherent bottom band most 
pronounced at k ~ (0, 0) the structure evolves with increasing 
either Ch or T into a single waterfall-like band which persists 
up to very high T ~ t or in the overdoped regime. 

The origin of the anomalous dispersion is according to our 
analysis entirely due to presence of strong correlations, as in- 
corporated already in the incoherent hole motion in a corre- 
lated insulator, as given within the Brinkman-Rice scenario 
Il24ll . Such a conclusion offers also the explanation why the 
waterfall phenomenon persists up to high T > J and in a very 
broad range of hole doping Cfe. Our results also indicate that 
explanations in terms of specific low-T features as the AFM 
long range order IU2I1 or AFM fluctuations Ilm[l7ll . quantum 
critical point ffl might be too narrow. In addition, similar 
waterfall behavior can be observed also in ID t-J model at 
high T > J H25I1 . However, instead of an incoherent bottom 
band a coherent holon branch emerges with reducing T < J. 



This is different from the present case in 2D. 

There are also some predictions relevant for ARPES ex- 
periments emerging from our analysis. In particular, ARPES 
spectra of electron-doped cuprates should correspond to uj > 
spectra of hole-doped cuprates (although with opposite t' 
and t") as already commented in 02611 . From the large asym- 
metry in uj as seen in presented results it follows that one 
cannot expect the HEK and waterfall in ARPES results of 
electron-doped cuprates. Further, our results predict an evo- 
lution of the anomalous dispersion with increasing T and c/j, 
nevertheless the waterfall features should persist up to very 
high T > J as well as in the overdoped regime. 
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